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Backscattering  From  Anisotropic  Random  Media 


1.  INTRODUCTION 

In  this  report  we  formulate  a  method  for  determining  the  backscattering  cross 
section  for  electromagnetic  (as  well  as  a  acoustic)  waves  from  slowly-varying 
anisotropic  fluctuations  in  a  random  medium.  In  the  electromagnetic  case  random 
media  having  these  properties  are  found  in  the  auroral  ionosphere  in  the  form  of 
striations,  that  is,  turbulent  plasma  that  has  become  elongated  into  thin  columns 
aligned  along  the  direction  of  the  Earth's  magnetic  field.  Similar  structures  have 
also  been  induced  artificially  by  barium  (Ba)  releases  in  the  ionosphere  above 
150  km  ’  and  by  high-intensity  ionospheric  heaters,  ’  The  method  developed 
here  is  based  on  the  cumulative  forward-scatter  single-backscatter  approximation 

5 

of  De  Wolf,  extended  to  include  anisotropic  media  as  well  as  backscattering  into 
(and  from)  arbitrary  directions. 

De  Wolf  employed  a  perturbation  expansion  of  the  Green's  function  for  the 
electromagnetic  wave  from  which  the  scatterecWield  and  cross  section  can  be 
determined  by  an  infinite  summation  of  selected  (bubble)  Feynman  diagrams,  which 
arise  upon  averaging  over  the  fluctuation  spectrum.  The  diagram  technique  for 
averaged  quantities  is  analogous  to  that  employed  in  quantum  field  theory  for 

(Received  for  publication  8  December  1982) 

Because  of  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  r'  >• 
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vacuum  expectation  values  with  the  correlation  of  two  fluctuations  replacing  the 
boson  propagator.  In  the  method  developed  here,  we  bypass  the  diagrammatic 
perturbation  approach  and  employ  an  operator  representation  for  the  full  Green's 
function  in  an  arbitrary  external  field.  In  this  approach  the  fluctuating  field  ap¬ 
pears  in  the  exponent  and  the  averaging  of  the  product  of  such  exponents  to  form 

the  expression  for  the  scattered  power  is  reexpressed  in  terms  of  correlation 

7 

functions  by  utilizing  Kubo's  cumulant  expansion  method. 

Actual  calculation  of  the  cross  section  for  a  given  turbulence  spectrum  re¬ 
quires  the  Fourier  transform  of  an  expression  involving  path  dependent  projections 
of  correlation  functions.  To  effect  this  Fourier  transform  we  utilize  a  technique, 

g 

similar  to  that  used  previously  in  the  calculation  of  multiphonon  absorption, 
which  involves  expanding  the  expression  to  be  Fourier  transformed  in  a  functional 
Taylor  series  in  powers  of  the  projected  correlation  functions.  The  subsequent 
Fourier  transform  of  each  term  can  then  be  carried  out  in  a  straightforward  way 
as  a  multiple  convolution  of  the  projected  correlation  functions. 

Calculation  of  the  backscattering  cross  section  has  been  carried  out  for  the 
case  of  a  Gaussian  correlation  function,  where  the  integrals  involved  may  be 
evaluated  analytically.  The  results  are  displayed  as  graphs  of  the  scattering 
cross  section  vs.  scattering  angle  for  various  choices  of  scale  lengths  and  thick¬ 
ness  for  the  random  media,  and  angles  of  incidence  and  frequency  for  the  incoming 
wave.  A  method  is  also  outlined  for  numerically  evaluating  the  cross  section  for 
the  case  of  the  more  physical  Kolmogoroff  spectrum,  which  essentially  involves 
representing  the  power  law  correlation  function  as  a  sum  of  Gaussians. 


2.  MULTISCATTERING  FROM  A  RANDOM  ANISOTROPIC  MEDIUM 

In  order  to  calculate  the  cross  section  for  backscattering  from  slowly-varying 
random  fluctuations  in  a  dielectric  medium,  we  will  initially  assume  that  depolari¬ 
zation  effects  are  small  and  that  propagation  through  the  medium  can  be  described 
by  the  scalar  wave  equation.  Denoting  the  propagating  E  field  by 

f(r,t)  =  E(?)  eiwt-£'r  (1) 


6.  Fradkin,  E.S.  (1966)  Application  of  functional  methods  in  quantum  field  theory 

and  quantum  statistics  (II),  Nucl.  Phys.  76:588. 

7.  Kubo,  R.  (1962)  Generalized  cumulant  expansion  method,  J.  Phys.  Soc.  Japan 

J£:1100. 

8.  Bendow,  B.  ,  Yukon,  S.P.,  and  Ying,  S,  C.  (1974)  Theory  of  multiphonon 

absorption  due  to  nonlinear  electric  moments  in  crystals,  Phys.  Rev.  B. 
10:2286. 
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the  wave  equation  for  a  given  transverse  component  E  may  be  written  as 

V2  E(r)  +  k2  e(r)  E(r)  =  0  .  (2) 

The  dielectric  permittivity  e(r)is  taken  to  be  equal  to  (e)  outside  the  scattering 
volume  and  equal  to  e(7)  =  (  e)  +  <5e(r)  inside.  The  scattering  volume,  as  indicated 
in  Figure  1,  has  been  taken  to  be  a  slab  of  infinite  extent  in  the  transverse  y-z 
plane  and  of  depth  L  in  the  direction  of  forward  propagation,  which  is  taken  to  be 
the  +x  direction. 


Figure  1.  Multiple  Small  Angle. 

Single  large  angle  scattering 
event  (a)  and  double  large  angle 
scattering  event  (b) 

In  the  following  we  neglect  the  contribution  of  reflections  and  backscattering 
from  the  interfaces  and  consider  only  backscattering  from  random  fluctuations 
within  the  slab  volume.  We  will  be  interested  in  calculating  scattering  from  irreg¬ 
ularities  that  tend  to  be  aligned  along  a  given  direction  such  as  a  background  mag¬ 
netic  field  B  (as  indicated  in  Figure  1).  We  will  continue,  however,  to  assume 
that  the  fluctuations  can  be  described  by  a  scalar,  albeit  one  that  is  in  general 
anisotropic. 


The  goal  of  the  present  calculation  is  to  determine  the  pou  t  r  scattered  out  of 
an  incoming  plane  wave  beam  of  unit  intensity  and  arbitrary  direction  of  incidence 
into  the  scattering  direction  K^.,  which  is  also  arbitrary.  The  definition  of  the 
angles  determining  IT  and  is  indicated  in  Figure  2.  Follow  ing  Ue  Wolf'’  we 
calculate  the  baekscattering  cross  section  on  the  assumption  that  the  scattering 
processes  that  dominate  are  those  composed  of  a  single  large-angle  hard  scatter¬ 
ing  event  followed  (and  preceded)  by  multiple  small-angle  forward  scatterings  ns 
indicated  in  Figure  1(a).  A  calculation  of  the  ratio  of  power  scattered  outside  a 
cone  of  apex  angle  2\c  to  that  inside  yields,  3  in  the  Born  approximation, 

I(\  -  \c>  "  (2kQLo  sin  (\c/2))"',/3  ,  Cl) 

where  kQ  =  u/c,  and  L  is  the  macroscale  of  an  assumed  3-IJ  Kolmogoroff  spectrum 
of  turbulence.  This  result  indicates  that  as  increases  beyond  the  small  angle 
(kQ l-*0 )  ,  the  portion  of  the  incident  beam  scattered  through  large  angles  becomes 

very  small.  A  similar  calculation  for  the  two  hard  scatter  processes  of  Figure  1(b) 
yields 

U\  >  \C)/I(\  c  \t.)  ~  (2kQ[,o  sin  (x c /  —  ))  6'  'J  (-4) 

with  even  smaller  ratios  for  high  order  processes.  Thus,  for  wavelengths  such 
that  kQLQ  >  1,  only  the  single  hard  scatter-multiple  small  angle  scattering  proc¬ 
esses  similar  to  Figure  1(a)  need  be  considered. 

We  define  the  scattered  wave  as  the  difference  betw  een  the  field  at  a  given 
point  r  with  fluctuations  present,  minus  the  field  in  their  absence.  Using  the 
Green's  function  G  for  the  perturbed  medium  that  satisfies 


V3  +  ^  ((  t)  +  de  (  r ) ) 


G(r  -  r  ')  =  6 ( r  -  r  ') 


and  for  the  unperturbed  medium 


■>  2  1 

f“  +  ~  <f  |  G°(r  -  r  ')  -  Mr  -  r  ') 


(a) 


(I  ) 


we  may  express  the  scattered  wave  as 


>  I  .  (  r) 


/ 


d  r 1  I  G  (  r  -  r  ’ ' 


G°(  r  - 


t  ’)!  Hr') 


(7) 


Figure  2.  Geometry  of  Initial  am!  Final 
Scattering  Angies 
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For  the  single  large  angle -multiple  small  angle  scattering  processes  of  Fig¬ 
ure  1(a),  the  series  of  Eq.  (9)  may  be  rearranged  to  give 


6E  (R) 
s 


■/ 


Gf(R,  r  ') 


fiefr ')  G.(r',R  )  dr’ 


(10) 


where  G.,  G^.  are  the  full  Green's  functions  for  propagation  along  K.  and  in  the 
multiple  small-angle  scatter  approximation. 

If  fie  is  assumed  to  be  a  normally  distributed  random  field  with  zero  mean, 
averaging  Eq.  (8)  leads  to  the  Dyson  equation  for  the  averaged  Green's  function  (G) 

(G)  =  G°  +  G°Q  (G)  (11) 

where  Q  is  the  mass  operator.  In  terms  of  Feynman  diagram's,  Q  is  represented 
by  the  sum  of  all  strongly-connected  diagrams® 


Q  = 


-J-X 


(12) 


where 

- - -  =  (6e(r)  6e(  r  '))  and  - - -  =  G°(r,  r  ')  . 

r  r'  r  r' 

2  2  2 

The  wave  number  shift,  6k  =  k  -  k  ,  due  to  the  random  fluctuations  is  then 

o 

given  by  the  Fourier  transform  of  the  mass  operator  Q  as 


9.  Tatarskii,  V.I.  (1971)  The  effect  of  the  turbulent  atmosphere  on  wave  prop¬ 
agation,  Israel  Program  for  Scientific  Translations,  Jerusalem. 


2 


(13) 


6  k2  =  /e*7 


Q( 7)  dr 


To  lowest  order  in  oe  only  the  first  diagram  in  Eq.  (12)  is  retained.  Assuming  a 
simple  exponential  for  the  correlation  function 


(6e(r  )  6e(r '))  =  e  7  7  ^ 


(14) 


leads  to  the  dispersion  relation" 


9  9 

k"  -  k" 
o 


x  2,4,2 
6  e  k  £ 
o 


,9  9  9  . 

Ik-  f"  +  (1  -  ikQ  i  )“  1 


To  lowest  order  in  6e,  for  kQf  S  1  this  gives 


(15) 


Re  (5  k) 


-k  6e2 
o 

8 


Im(6k) 


a  2,2, 
6e  k  f 
o 


(16) 


thus  there  will  be  both  a  wave  number  shift  as  well  as  attenuation  due  to  multiple 

scattering  in  the  random  medium.  As  an  example,  a  30-MHz  signal  traversing  a 

random  medium  with  1  percent  permittivity  fluctuations  having  a  100-m  scale 

6k  -3 

length  would  experience  a  wave  number  shift  — =  1.25  X  10  and  an  attenuation 
of  a  =  1.0/  km . 

To  calculate  the  scatter  cross  section  in  the  Born  approximation,  we  approx¬ 
imate  G.,  Gj  by  G°,  where 


G°(7,  ?') 


1 

T7 


(17) 


The  expression  for  the  scattered  field  at  R  for  a  delta  function  source  at  Rq  given 
by  Eq.  (10)  may  then  be  written  as 


6E°(R) 

s 


G°(R,  7')(-k2 


6eCr'))  G0(7,.Ro) 


(18) 
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In  the  plane  wave  limit  (H  ,H  -  -o)  we  may  approximate 


where  m  is  the  angle  between  the  polarization  direction  of  the  incident  wave  E  and 

-*■  •>  ° 
the  scattering  direction  K^.  In  what  follows  we  shall  leave  off  the  sin"  (u.)  polari¬ 
zation  dependence  and  consider  only  the  scalar  quantities  i_  ^  and  o_.  _  . 


Thus,  for  the  Born  approximation  we  have 


Ki'Kf 


K.,  K, 
1  f 


..2  r  i(K  -K.)-  r  1 

fK  £  =wjdr'e  6e(r,) 

i’  f 


Ki*Kf 


=  /dr'/dr' 


i(K.-K. )(?'-?") 

e  (6e(r')  6e*(r  ")) 


Changing  to  new  variables  Ar  =  r  '  -  r  ",  r  =  - g -  ,  and  AK  =  K.  -  K^.,  this  may 

be  written  as 


KiiK[ 


K4  e2  V 


4-  (AK)  , 


where  4>  (K)  is  the  Fourier  transform  of  the  correlation  function  for  a  homogeneous 


medium 


-i  (6 e2(?  -?'))  =  -ir  (6e(r  )  6e!':(r  ')N 


We  note  that  if  (<5e  tAr))  is  Gaussian,  then  $  (K)  is  also  Gaussian.  We  shall 
assume  that  4>  can  be  factored  as 


<MKx,Ky,Kz>  =  d(Kz)^(Kx,Ky)  . 


(28) 


The  two  functional  forms  that  are  of  greatest  interest  are  the  Gaussian 

_  _  -  _  -K2  L2  / 4I  T  ,  -K2L2/4 

®  M  (K)  =  d(K  >j9~(K  )  =  S*  L  e  z  11  »L2  e  1  -1 

47  2  Z7  1  L  '!  J  L  i 


^2  ^ 

*^(K)  =  d(K  )i(K)  =  e  z  11 

z  rt  1  II 


1  +  K2  L2 

L  1  ij 


M+l  ’ 


where  m  =  5/6  corresponds  to  the  Kolmogoroff  spectrum  for  turbulence.  In  the 
following  we  shall  assume  that  L ^  >  L^. 

To  include  multiscattering  effects  in  the  simplest  way  we  utilize  the  eikonal 

Q 

approximation  for  G  given  by 


iy  /  6e (r +s'K)ds' 


G(7,f ')  =  G°(r,r ')  e 


where  the  integral  is  to  be  taken  over  the  effective  path  length  in  the  medium.  If 
we  restrict  the  possible  scattering  angles  to  those  outside  a  small  cone  of  half 
angle  ec  -  around  the  direct  backscatter  direction  (Kf  =  -K.)  then  the 

averages  over  products  of  G.  and  Gf  are  essentially  disjoint  and  may  be  factored. 
Using  the  eikonal  approximation  and  grouping  terms  on  each  path  together,  we  get 


d?"/d?'  eiAg-(?'-?") 
K.,Kf  4v  J  J 


f  iy  f  <5  e  ( r '+cjK  .  )do  -iy  /  6e’|!(r  "+a'K.)da' 

..  /  -ad  1  -oc  1 

X  \  e  e 


X  (6e( r  ')  6e* ( r  ")) (e 


iy  J  6e*(1r '+aKj)da  +iy  J  6e (r  "+a'Kj)da’ 


16 


,.4  r,~  /*..-»  i-i K -  Ar 

a _  _  =  K_  I  dr  dii-  e 

K.,  Kj.  4ir  J  J 


K~S.  x 


j  ({6f"  (Ar  +  1vt)n  -{6c ~  (K^t)n  }  dr 


>  e 


K'sf  *  - 

— ^ —  /  { ( 6 £ “  (Ar  +  Kp))  -(6e“(K^t)n  }Ut 

X  (6e“(Ar))  e  00 


(37) 


We  now  reexpress  the  path  dependent  terms  in  the  exponent  in  terms  of  their 
Fourier  transforms  as 
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(6  e"  (Ar  +  K.t))  dr 


iq.  (Ar  r  K.t)  ^  9 
e  ',(if"(q)N 


kV  f  df 

'*  J  (2,r 


t.iq-  4? 


q>  *  (q) 


(38) 


i'h(  delta  function  6(f..  •  q )  a.-ts  os  ;■  pro..*-.  ti«-n  op-,  vjioi,  -liowi  'g  c  niv  those 
.  t.  of  q  that  lie  in  tin  plane  pcrpe:.'h<  ulat  to  K  j’.  ‘ht  origin  "o  -ontribui.e  ‘o 
*  hi;  i  o..ri*.r  integral .  ‘A  «•  is*  thus  1  to  defi.  ••  ‘he  (pa'h  !<  pen  uenl'i  ;  ■  n’t  •  f «  i 
"it.  .  "laiion  tune:.  > ■: .  I  -  (Ar)  ihro.teh  ■  r<  In'  .■ 


Performing  the  dr  integration  yields 


* 


b 


and  the  iterated  convolution  through 


Dn+1(K)  =  f  Dn(S  -p)D(p) 

J  (2*r 


(51) 


The  double  Taylor  series  expansion  of  the  functional  (ef  -  l]/f  about  D.(r), 
Df(r)  =  0  thus  yields 


oc  n 


n=0  j=0 


ef-  1 


Dn~j(Ar)  D?.<Ar)  .  (52) 

|f=-(0.+3f)L  1  f 


n  f 

The  derivative  S^de1  -  lJ/f)|j__x  may  be  expressed  in  terms  of  the  incomplete 
gamma  function  yin.x),  as  (see  Appendix  B) 


c.n/r„f  <h/(l  _  yin  +  1),  x) 

9fUe  l])/f)|f=.x - T+T - 


(53) 


The  final  expression  for  the  scattering  cross  section  for  a  unit  volume  is  then 
given  by 


Eln]  K4  f  iAK- 
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°i'l'f  n=0  Ki,Kf 
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X  J  pJ.  - i— -i -  •  Dn  J(Ar>  DJr(Ar) 

Hf  .  /o  .  ,  m+l  l  f 


Ujh  +  4lf)n+1 


(54) 


where  the  n  =  0  term  of  the  infinite  sum  has  been  separated  out  as  the  first  term  in 
the  large  brackets,  and  is  essentially  the  Born  term  modified  by  the  attenuation 
along  the  forward  and  backward  paths. 

The  n'''1  convolution  of  the  projected  correlation  function  defined  in  i  qs.  (38) 
and  (3d)  is  calculated  in  the  most  straightforward  manner  by  using  the  delta  func¬ 
tion  in  K q .  (33)  to  integrate  over  the  q  coordinate  paraih  1  to  K.  and  then  dcvelopim. 


the  iterated  (2-lJl  convolution  series  as  defun  <1  in  l.q.  (>]'  a  tin  n  ,  ,  r  a 

ular  to  K.,  Hi  is  approach  runs  into  ultie:-  due  t,,  tin.  :u.  t  that  iuu  : 

functions  of  Kqs.  (29)  and  (do)  factor  into  functions  defin<  <1  a  to  nr  /  ;.i>.  :  ’in 
plant-  ratiier  than  K  (or  K^,)  and  tile  plane  perpendicular  to  K.  (or  id,.),  u.  a  -  -tUo: 
the  general  term  in  the  sum  of  ICq.  (54)  requires  the  <  onvolution  of  L>‘‘  J  and  1/ 
which  are  defined  on  planes  of  differing  orientation.  The  simplest  way  to  avoid 
both  of  these  difficulties  is  to  use  the  projection  delta  functions  in  Tq.  (otij  lo 
integrate  over  the  q^  coordinate,  projecting  both  D1.1  •*  and  L)j.  onto  the  q.  -  q 
plane. 

To  actually  calculate  a  given  convolution  product  it  is  best  to  form  the  (n  - 
convolution  of  L).  in  a  coordinate  system  p  ,  p  ,  where  p  is  directed  along  K; 

i  .  j.  fl  '  t 

and  where  is  the  projection  of  K.onto  the  p  -  p  plane  as  shown  in  Figure  ;i 

j.  1  ~  j  y 

In  a  similar  manner,  the  j  convolution  of  Dj.(q)  can  be  formed  in  a  coordinate 
system  where  q  ^  is  directed  along  Kf  . 


Figure  3.  Decomposition  of  Scattering  Wave  Vector  Into  Com¬ 
ponents  Parallel  and  Perpendicular  to  Magnetic  Field  Direction 


One  can  show  that  the  convolution  of  two  K  projected  functions  is  another  K 
projected  function  (see  Appendix  C).  Thus,  for  the  Gaussian 


(55) 


•)0 


we  get  for  the  projected  correlation  function  in  the  px  -  p  plane  a  hybrid  given  by 


D^‘j(p)  =  D”'j(p  )[2JT6(K.  ■  p)  K .  cos  (\ . ) | 
i  i  1  *  l  i 


(56) 


with 


’rLlI'i  -PvL"/4(n-j)  -p“  L^/  4  (n-j) 

TFT^-jy-  e 


<57 ) 


where 


Li  =  LX  +  Ll!  tan“’  (;M)  * 


(58) 


A  general  term  in  the  sum  of  Eq.  (54)  with  (n  -  j)  >  0,  j  >  0  will  require  the 
Fourier  transform 


J dAr  e^K‘Ar  <6e2(A~))  D^Ar)  dJ(aT)  , 

which  we  reexpress  as  a  double  convolution 


(59) 


r  _  i(K.-Kf)-Ar  r 
J  dAr  e  1  f  / 
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Ar 
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dpdqdf  ,  .  ~_-i  —  -  —  - 

~ '  fi -  <<5f“(f  )'  D.  J(p)  DJ(q)  6(f  +  p  +  q  +  K.  -  Kf)  . 


(60) 


(-hanging  variables  to  p*  ’  =  p  +  K.,  q'  +  -q  +  Kf  and  dropping  primes  this  becomes 
after  integration  over  d  ? 


fi '  V 


:  6  (  (p  -  K.)  •  K.  )  .  6  (  ( q  -  Kf)  •  Kf  )  (2trK.J(2trKf ,)  . 


This  can  now  be  evaluated  by  integrating  the  two  deita  functions  over  p  and  q^ 
and  projecting^.  along  the  p  axes  as 


P1  ‘  Iq.J  =  cos  ‘  '  ^9  sin  _  Oj) 


p9  •  [ "q±J  =  q  sin  (of  -  +  q2  cos  (of  -  Co) 


to  form  the  argument  for  (6e“(q  -  p)) .  These  operations  are  carried  out  in 
Appendix  D  for  the  Gaussian  correlation  function  of  Eq.  (2  9). 

The  terms  in  the  sum  of  Eq.  (54)  with  (n  -  j)  =  0,  j  >  0,  and  (n  -  j)  >  0,  j  =  0 
can  be  evaluated  in  a  similar  manner.  Results  for  the  Gaussian  correlation  func¬ 
tion  are  developed  in  Appendix  E.  The  expression  for  the  total  cross  section  for 
the  Gaussian  correlation  function,  using  the  results  of  Appendices  L>  and  E,  thus 
becomes  a  simple  double  sum  of  exponentials.  This  expression  has  been  evaluated 
for  various  values  of  the  medium  parameters  L  ,  L  ,  slab  thickness  L,  wave 
number  K,  and  angle  of  incidence  (and  is  plotted  in  Figures  fia  through  lie). 

The  behavior  of  the  cross  section  for  the  Gaussian  correlation  function  is  discussed 
in  the  next  section. 

The  evaluation  of  the  cross  section  given  by  Eq.  (54)  for  the  Kolmogoroff 
spectrum  of  Eq.  (an)  is  more  difficult  due  to  the  fact  that  integration  of  the  projec¬ 
tion  delta  function  of  Eq.  (38)  over  p  leads  to  an  expression  for  lJl*fj  which,  as 

2  „  ■> 

in  the  Gaussian  case,  is  no  longer  a  function  of  p  7  but  of  p  , ,  and  p  separately 
and  which,  moreover,  does  not  factor  as  does  the  Gaussian.  The  2-1)  convolutions 
for  bt^C|  can  be  carried  out  numerically  in  a  straightforward  way  with  no  difficulty 
other  than  that  of  the  computation  time  and  core  memory  involved,  since  the  scries 
of  Eq.  (54)  converges  rapidly  (N  =  5  is  sufficient  for  the  Gaussian)  this  is  not  a 
severe  problem  unless  a  large  number  of  points  (for  example,  as  represented  bv 
the  curves  of  Figures  5a  through  lie),  are  involved.  In  order  to  compute  o  for 


the  Kolmogoroff  spectrum  in  a  more  economical  way,  we  represent  the  correlation 
function 
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d + 1 L! 
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as  a  sum  of  Gaussians  by  employing  the  integral  representation  for  the  gamma 
function 
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(64) 


where  w.  and  x.  are  the  Gauss -Laguerre  weights  and  zeros.  This  allows  the  func- 

v  t  t  v 

tion  D[*J  to  be  approximated  by  a  sum  of  separable  Gaussian  Dl^Cl  as  given  in 
Eq.  (57).  A  way  to  avoid  retaining  ail  of  the  (M)n  terms  that  would  appear  in  the 
n  convolution  of  Di^J  is  to  approximate  each  n  order  convolution  Ul^l  by  an 
appropriate  new  sum  of  Gaussians.  A  procedure  for  doing  this  is  given  in 
Appendix  F. 

Although  we  have  utilized  the  eikonal  approximation  [Eq.  (5)1  to  the  full  Green's 
function  in  the  calculation  thus  far  it  is  possible  to  give  a  completely  general  for¬ 
mulation  for  the  cross  section  by  using,  in  conjunction  with  the  cumulant  and  con¬ 
volution  expansions  of  the  previous  section,  the  operator  representation  for  the 

n 

full  Green's  function  in  an  arbitrary  external  field  developed  bv  Fradkin  and 


ti 


Milekhin.  ^  In  this  method  the  full  Green's  function,  starting  from  the  Fourier 
transform  of  Eq.  (5)  is  given  bv 

G(r,p)  =  -  [(v-  ip)2  -  k2  e (r)j 


OC 

■*/ 


is  (V- ipE-k^ef  r)j 


ds 


=  ije 
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is(-p2+isV2+2  p- V- isk  e(r) 


ds 


(66) 


Milekhin  has  shown  that  the  principal  approximation  that  is  necessary  in  deriving 
the  eikonal  limit  corresponds  to  dropping  the  term  [is  v“i  in  the  exponent  of 
Eq.  (66).  Operating  with  e2sp’  ^  on  e(7),  shifts  the  argument  from  r  to  7  +  2ps 
and  operator  ordering  leads  to  the  eikonal  form  given  in  Eq.  (6), 

Calculation  of  the  scattering  cross  section  to  include  the  effects  of  terms  in 
the  Green’s  function  beyond  the  eikonal  approximation  is  most  straightforward  if 
Fradkin's  representation  using  functional  derivatives  is  employed  [Reference  6, 
Eq.  (5.  17)] .  We  shall,  however,  leave  the  calculation  of  the  cross  section  includ¬ 
ing  these  higher  order  contributions  to  a  future  publication. 


3.  CROSS  SECTION  FOR  GAUSSIAN  FLUCTUATIONS 


G:v- 


The  behavior  of  the  cross  section  can  most  easily  be  studied  for  the  case 
where  the  correlation  function  (6e2( r  -  ?'))  is  given  by  the  factored  Gaussian  of 
Eq.  (29).  Looking  first  at  the  n  =  0  (modified  Born)  term  in  Eq.  (54)  we  have 


r- 
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1  or  L  »  1.  this  tern;  'A  ill  be  .small  unless  i..  K,  1  in  Iku  ::  t<  ..:i  ..  hi  thus 
peak  when  \f  ~  \.  leading  to  essentially  specular  seam -ring.  Mu'uscatu-!  mg  effects 
are  contained  in  the  term  in  brackets  and  appear  as  a  path  dependent  attenuation 
for  traversing  the  slab,  f  or  small  1.  the  exponent  can  be  expanded  and  the 
bracketed  term  equals  ~  1,  which  is  expected  as  multiscatter ing  effects  will  not 
contribute  unless  the  path  length  is  longer  than  the  mean  free  path  between  small 
angle  scatterings  (<£  l/2«).  Due  to  the  anisotropy  of  the  medium,  paths  direc¬ 

ted  along  the  x  axis  I  X  ~  ^  ]  will  be  attenuated  less  per  unit  length  than  those 
directed  more  toward  the  Z  axis  (x~-  0)  as  can  be  seen  from  the  expression  foru 
given  in  Eqs.  (41)  and  (42).  Recalling  that  gh  contains  the  geometric  factor 
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the  product 


,j.  L  =  -la  ■  !.  •  T(K.) 
■i  i 


->) 


is  seen  to  represent  the  path  length  from  the  front  of  the  slab  to  the  back  of  the 
slab  along  the  rav  K;  times  the  attenuation  lor  that  rav  direction.  Thus  for  shallow 
angle  rays,  In  addition  to  the  larger  attenuation  as  a  function  of  the  polar  angle  x> 
there  will  be  an  added  reduction  in  the  scattered  power  due  simply  to  the  longer 
path  lengths  that  are  required  in  traversing  slab  at  shallow  angles. 

In  order  to  determine  the  behavior  of  the  n  >  1  terms  in  Eq.  (54)  in  an  approx¬ 
imate  way,  we  shall  assume  that  L  »  E  so  that  the  scattering  will  be  approxi¬ 
mately  spectral  and  set  p.  Taking  (T  L  >  1  we  approximate  -y  (n  +  1,  2p.E) 

by  I’ (n  +  1)  n!  i'he  terms  with  n  >  1  are  then 
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Using  Eq.  (57)  for  the  n^1  order  convolution,  we  then  have  very  roughly  an  approx  - 

[n  I 

imation  lor  the  behavior  of  a  given  by 
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The  contributions  to  a  from  the  n1^1  order  multiscattering  terms  are  thus  re- 
K.K„ 

if  _  1 

duced  in  amplitude  by  a  factory  (n  +  1)  and  are  broadened  by  a  change  in  the 
effective  half  width  from  l/L  ^  to  n  +  l/ L(|.  This  behavior  is  illustrated  in  the 
sketch  of  Figure  4. 


Figure  4.  Sketch  of  Expected  Behavior  for  o  (Born)  and  nth 
Order  Scattering  Cross  Section  Contributions 


f  n  I 

Since  we  have  indicated  that  each  term  al_  decreases  in  amplitude  as 

K..Kf 


N 


(n  +  1) 


the  series 


In] 


would  appear  to  diverge  as  lim^^^  ln(N).  However, 


,  K..K, 
n=l  i  f 

if  the  approximation  y(n  +l,x)  ~  T(n)  is  not  made,  convergence  of  the  series  can 
be  proved  by  using  the  integral  representation  for  y{n  +  l,x) 


y(n  +  1,  x) 
n! (n  +  l) 


n  =  0 


00 


L 

n=0 


/ 


-t  .n 
e  t  at 

(n  +  l).r 


o  n  =  0 


1 


n+1 


(n  +  1).' 


dt 


x  ♦  r  00 

r  - 1  x — \  j.tn 


■/4  E5 

o  _  m=0 


-E.(-x)  +  C  +  in(x) 


E1(x)  +  C  +  l  n(x)  , 


where  E.(x),  E^fx)  are  the  exponential  integral  functions  and  C  is  Euler's  constant. 
For  the  case  of  an  isotropic  medium  L  =  and  a _  will  no  longer  depend 


K^Kf 


on  and  individually,  but  only  on  their  difference  Alt  =  ^  -  K..  The  momen¬ 
tum  transfer  AK  may  be  written  as 


AK2  =  (K.  -  Kf)2  =  K2(2  -  2  cos  ( 6 ))  =  K2  2  sin  ( 


where  0,  the  angle  between  K.  and  Kf,  is  defined  by 

cos  (0)  =  FT  •  Kf  =  sin  (Xj)  sin  (x^)  cos  (< f>.  -  $f)  +  cos  (Xj)  cos  (x^) 

where  the  angles  x,  <t>  are  defined  in  Figure  2.  If  we  take  =  0,  <t>^  =  r  then 
0  =  Xj  +  and  we  thus  expect 


Ki(Kf 


-AK2*2/4(n+l) 
(n  +  1) 


to  decrease  (for  Xj(>  0)  fixed)  as  \j.  increases,  as  sketched  in  Figure  5  for  the  case 
L  =  0.  For  a  slab  of  small  but  finite  thickness,  the  path  length  in  the  Xj  =  0 
direction  becomes  infinite,  cutting  off  scattering  in  this  direction,  whereas  for  a 
thick  slab  the  scattering  will  be  reduced  for  all  Xj  >  0  as  well.  This  behavior  is 
sketched  in  Figure  5  where  <  0  means  =  jxjl.  ^  =  0,  that  is,  scattering  into 
the  forward  direction. 


* 

ft 


as  functions  of  x*  for  x{  fixed  (with  <t>.  =  0,  <t>f  =  jr,  e  =  0.  03)  The  starred  curve 
fol  -1  1  11 

represents  a^J  ,  the  highest  solid  curve  a_„  ,  and  their  difference  is  the 

remaining  solid  curve  a 


In  Figures  6a  through  6e  we  have  plotted  the  cross  section  for  the  case  of  a 
thin  slab  (L  =  l. )  for  which  the  total  cross  section  essentially  equals  the  Born 
cross  section  with  multiscattering  effects  being  mainly  due  to  the  attenuation  of 


the  random  medium,  as  discussed  as  the  beginning  of  Section  3.  The  anisotropy 
L 

aspect  ratio  p  =  ^ 11  varys  from  p  =  2,  in  Figure  6a  which  appears  almost  iso¬ 
tropic,  to  p  =  32  in  Figure  6e.  The  expected  behavior  for  p  =  1  was  discussed 


previously  and  sketched  in  Figure  5.  For  r)  =  2  the  situation  is  seen  to  be  quite 


similar.  There  is  a  small  contribution  due  to  a 


M 


Ki,Kf 


which  falls  off  exponentially 


as  increases.  As  17  is  increased  to  tj  =  4  the  Born  term  develops  two  peaks. 
From  Eq.  (67)  we  see  that  the  Born  term  will  be  given  by  a  product  of  Gaussians 
(for  the  case  L  =  1  the  attenuation  factor  is  ~  1) 


!°i  ,  1 


_  =  f  e2  K“ 
Ki-Kf  4 


2  -4K2iA,2/4  -ai?L2/4 

L!  S  *  22  .  22 


with 


AK2  =  K2  (cos  (Xj)  -  cos  (xJ)2 


AK2  =  K2  (sin  (Xj)  +  sin  (Xj))2 


and 


cos  (Xj>  =  sin  (\j)  =  ,  xi  =  45° 


->9 

The  Gaussian  with  argument  proportional  to  AK“  peaks  at  xf  =  0,  while  trat  with  ar 

2  2  ^  T 

gument  proportional  to  AK“  n  peaks  at  Xj.  =  X-  Their  product  will  thus  have  a  c'ir 
imum  somewhere  between  =  0  and  x^  =  Xj  (unless  p  =  1)  and  this  is  reflected  by 

the  dip  seen  in  Figure  fib  at  ~  11°.  As  p  is  increased  to  p  =  8  in  Figure  fic,  the 

2  ■>  •>  1 

term  exp  (*E1  K“  p“ /4 )  dominates  and  the  scattering  cross  section  becomes  peaked 
around  Xj-  =  Xj* 

Increasing  p  to  p  =  lfi  and  p  =  32  as  shown  in  Figures  fid  and  fie  results  in 
narrower  peaks  with  larger  maxima.  The  integrated  cross  section  remains 


approximately  constant  however,  since  the  correlation  function  is  normalized, 
that  is. 


vV  L , 


00 

/ 


-K2  L2  / 4 

dK  e  2  "  =1 

z 


and  for  very  narrow  peaks 


dK 


<z  =  d[K  cos(\f)|  =  -(sin(xf)  Kj  dyf  ~  const  •  d\f 


as  sin  (\j.)  is  essentially  constant  over  the  main  contribution  to  the  integrated  cross 
section. 

In  Figures  7a  through  7e  ail  parameters  remain  the  same  as  in  Figures  6a 
through  6e  except  that  the  random  medium  is  now  taken  to  be  a  thick  rather  than 
a  thin  slab  (L  =  10000  vs.  L  =  1).  Multiscattering  and  path  length  attenuation  will 
thus  play  a  larger  role.  Comparing  Figures  7c  through  7e  with  Figures  6c  through 
6e,  it  can  be  seen  that  the  change  in  the  modified  Born  contribution  _  is 


entirely  due  to  the  attenuation  factor 
-o 


-(/i.-t-pfJL 


Ki,Kf 


1  -  e  J/  [(|F  +  Mf)LJ  .  which  around 

the  peak  Xf  =  45",  equals  ~  0.0046.  The  multiscattering  contribution  for  the  thick 
slab  case  is  broader  than  the  Born  term  as  discussed  previously  and  is  of  com¬ 
parable  magnitude.  For  Figures  7a  and  7b  the  implication  of  a  thicker  slab  for 
the  Born  term  is  that  the  attenuation  factor,  instead  of  being  ~  1  for  all  angles 
will  now  be  angle  dependent,  being  much  smaller  at  shallow  angles  due  to  the  in¬ 
creased  path  length  in  the  medium  (as  well  as  the  effect  of  a ^  being  somewhat 
larger  for  shallow  angles).  This  has  the  effect  of  suppressing  the  small  angle 
peaks  of  Figures  6a  and  6b  resulting  in  the  broadened  and  reduced  peaks  of 
Figures  7a  and  7b. 


The  multiscatter  contributions  cr 


M 
K.,  Is.. 

i  1 


in  Figures  7a  and  7b  are  seen  to  be 


broader  but  roughly  similar  in  shape  to  the  modified  Born  terms.  They  are,  how¬ 
ever,  of  greater  magnitude  being  more  than  twice  the  Born  cross  section  for  q  =  2. 

In  Figures  8a  through  8e  has  been  reduced  from  x^  =  45°  to  x^  =  15°.  all 
other  parameters  being  the  same  as  those  in  Figures  7a  through  7e.  The  main 
difference  in  the  plots  for  large  rj,  q  =  8,  16,  32  is  that  the  peaks  for  xi  =  15°  ire 

much  broader  than  those  for  x;  =  45°  about  their  respective  maxima,  as  well  as 
2  1 

being  larger  (by  ~  10  ).  If  we  look  at  the  dominant  factor  in  the  Born  term 


b.(K.  )  =  e 
l  fz 


-(K.  -K.  )“  I ,  / 4 
iz  fz 


,  this  would  yield  a  curve  with  the  same  half  width  about 


any  value  for  K  if  plotted  as  a  function  of  K  .  As  f'gurt  -  >  a  Uw  m,ci:  1  1<  have 
been  plotted  as  a  function  of  y^.,  a  plot  of 

b  lK  COS  ( yd  I  vs  x.^. 

will  be  broader  for  curves  peaked  around  smaller  y^  since 
idK„  i  =  siniy.)  K  •  dyf 

that  is,  equal  increments  of  dy,.  yield  a  smaller  change  in  for  small  y,.(\j-  "  0°) 
than  they  do  for  large  \r(\..  '  up'1).  The  fact  that  the  values  of  a  . '  .  at  w  =  la” 

i;i'Kf 

in  figures  lia  through  tie  are  larger  than  in  figures  7a  through  7c  is  due  to  the 
second  factor  in  tin-  Horn  term  lfq.  f'V)| 

-AK2I  2  ■  !  -  (sir.  \.-sin  \,ri-' !  ~/4 

1  1  ! 
t  / 


-(>.  ■  •  To,  r  !.-  j.-  '4 


•!  .  >  -•  i 


r  f  /  m 


aiong  with  a  factor  of  -•  o.  12a  clue  to  tnc  larger  attenuation  for  the  shallow  angle 
paths  in  figures  fic  through  8e,  lead  to  the  observed  difierencc  in  these  peak 
maxima  which  amounts  to  a  factor  of  '  fi  with,  respect  to  the  peak  maxima  of 
Figures  7  c  through  7  c. 

for  x-  =  4 5°  the  peak  multiscattering  contribution  with  respect  to  the  Horn 


contribution  runs  from  about 


-  ii.  i  lor  r>  =  82  to 
l°i 


~  2. 3  for  r,  -  2 ,  while 


foi  y.  -  1  7  the  ratio  runs  from  m  •;  to  i- 

in  figures  Pa  through  he  tin-  •toss  so-  lion  for  rj  =  4,  y.  =  45°,  K  =  2  *•/ 1 .  (5  is 
plotted  for  I.  r  1.,  id,  Thi,  loon,  and  lOOOu.  The  double  peak  in  figure  Pa  was 
discussed  previously  for  f  igure  fib  (f  igure  f>b  =  figure  Pa).  As  1.  is  increased  to 
1  n  the  low  angle  peak  in  the  Horn  term  is  suppressed  whiic  the  peak  around  \ .  is 
lowered  slightly,  fince  the  only  1.  dependence  in  the  Horn  term  |f.q.  (f:7)|  is  in  the 
attenuation  factor,  the  low  angle  suppression  is  due  to  the  increase  in  g,.  at  low 

As  i.  is  increased  beyond  10,  the 


angles  (for  1=1  the  attenuation  factor  is  ~  !). 


attenuation  factor  behaves  roughly  as  l/(j i.  +  L,  which  is  reflected  in  the 

[o )  *  * 

behavior  of  a_ ^  in  Figures  9c  through  9e. 

k  . ,  i\ 

i  i 

In  Figures  10a  through  lOe  and  11a  through  lie  we  have  plotted  the  multiscat¬ 
ter  cross  sections  for  increasing  values  of  x^»  x^  =  5°,  15°,  45°,  75°,  85°;  for 
K  =  2 w/  1.6  (n  =  4,  L  =  10000)  in  Figures  10a  through  lOe,  and  for  K  =  2?r/0.8 
(r?  =  4,  L  =  10000)  in  Figures  11a  through  lie.  In  Figures  10a  and  10b  the  peak 
cross  sections  are  roughly  comparable  and  _  is  roughly  equal  to  ^  . 


Ki'Kf 


K K 
i  f 


For  =  45°,  7  5°,  and  85°  (Figures  10c  through  lOe)  is  still  roughlv 

K.,K. 

(ol  1  * 

comparable  to  cr_  _  however  the  peak  cross  sections  drop  rapidly  for  X-  =45° 


Ki'Kf  -AK^L^/4 

and  beyond,  due  to  the  factor  e  1  ,  which  becomes  exponentially  smaller 

around  the  peak  at  as  —90°.  Thus  a  doubling  of  K  from  K  =  2»/  1.  fi  to 
K  =  2  jr/o.  8  means  that  the  cross  section  around  the  peak  at  x^  =  X^  is  given  roughly 
by  ct(K  =  2n/ 0.  8)  ~  (ct(K  =  2?r/l.  6)]  which  accounts  for  the  orders  of  magnitude 
difference  between  Figures  10c  through  lOe  and  the  comparable  Figures  11c 
through  lie.  If  instead  of  the  Gaussian  correlation  function  we  had  employed  the 
Kolmogoroff  spectrum  [Eq.  (30)  withp  =  5/6],  this  would  lead  to  a  corresponding 
change  in  the  cross  section  around  the  peak  of  roughly 


that  is,  one  to  two  orders  of  magnitude  vs.  orders  of  magnitude  for  the 

Gaussian. 
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Figure  7d.  Plot  of  Modified  Horn,  Mult  [scattering,  and  Total 
Cross  Section  as  a  Function  of  x.-  for:  x,  =  15  ,  rj  =  16, 

L  =  1  uDilii,  >  =  1.6  1 


ANGLE 


Figure  7e.  Plot  of  Modified  Born,  Multiscattering,  and  Total 
Cross  Section  as  a  Function  of  vf  for:  =  4  5°,  rj  =  32, 

1.  =  10000,  x  =  1.0  1  1 
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3i 


ANGLE 


Figure  10c.  Plot  of  Modified  liorn,  Multiscatjering,  and  Total 
Cross  Section  as  a  Function  of  \,  for:  \.  =  4f>  ,  rj  =  4, 

I.  =  10000,  >  =  l.o  1 


Figure  lOe.  Plot  of  Modified  Born,  Multiscattering,  and  Total 
Cross  Section  as  a  Function  of  x-  for:  \.  =  85  ,  n  =  4, 

L  =  101)00,  A.  =  1.(5  1  1 
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Appendix  A 

Calculation  of  Effective  Path  Length 


The  path  length  AC  may  be  determined  in  terms  of  the  polar  and  azimuth 
angles  \,  <p  of  the  ray  direction  K,  and  the  distance  into  the  medium  x  by 

AC  =  (AB2  +  BC2)1/2  ,  (Al) 

where 

x  =  BD,  AB  =  x  tan  and  BC  =  x/sin  (\)  (A2) 


as  indicated  in  Figure  Al. 


A 


Figure  Al.  Geometry  for  Determining  Effective 
Path  Length  for  a  Given  Scattering  Wave  Vector 


Appendix  B 

Evaluation  of  3 ^  [e^  -  1]  /f  in  Terms  of  the  Incomplete 

Gamma  Function 


th  f 

The  n  derivative  of  (c  -  l)/f  with  respect  to  f  may  be  reexpressed  by  using 
the  Leibniz  theorem  as 


f  e1  -  1 

1 

zs 

1 

1 

L  f 

■  '  LTJ 

-  o 

.  i  . 

--  |  oVxr1)  +  (^xa"-1  efX3f’1) 

+  0<3n-2  <■'  ho1  r1)  *  . . .  +  e1 3na‘ 1 ) ;  -  o»(f'i)  .  (BJ 

1  :sinu 

",(l  ^  ,  1  (f"  1  1  -  >t"\  ••  -Ar1,  = 


OVM 


this  becomes 


i  i  +  (j>(-f_1)  +  (2 > 2 i  +  •  •  •  +  (")j:  (-f'1)3 

+  •  •  •  +  n!  (-f-1)" 


thus 


ef  -  1 


fx-rM  4  - 

i  j  n.  i 

j  =  0 


{f-l)n 


(B3) 


-eLY  _ L_  (.f-ijj  -1 

f  2-j  (n  -  j)!  (  f  ’ 


-f-V 


3  =  0 


Letting  k  =  n  -  j  this  can  be  reexpressed  as 


(B4) 


i-  E  irr<-t‘1,n"k  -  ^  ■  4  E  ilr  '-°k  - 


,k  (-f_1)n 


k=n 


k=0 


(-I)n  j  j  V'  (-f)k  .  A 

fn+l 


E-4# 


k=0 


•?  ' 


Evaluating  this  at  the  point  f  =  -x  yields 


[V-il 

.  -n! 

"  k 

*'x E vr-‘ 

f 

,  n+1 

1  f  =  -X  X 

k=0 

(B5) 


(Bfi ) 


The  series  expansion  for  the  incomplete  gamma  function  is  given  by 


y(n  +  1,  x)  =  n! 


(137  1 


7b 


Appendix  C 

Convolution  of  Projected  Correlation  Functions 


To  show  that  the  convolution  of  two  projected  correlation  functions  is  another 

vn 

projected  correlation  function  we  start  from  the  definition  of  O^Cp^)  in  Eq.  (5fi), 
where  IE(p  )  is  defined  in  the  plane  p  -  p  perpendicular  to  p  at  the  origin 
<PZ  ■  0). 

D.(p)  =  D.(pi)[2jr6(p  •  K)Kz] 


=  D.(pi)27r6 


(Cl) 


The  convolution  of  D.^  with  D.  to  yield  D?  is  defined  in  Eq.  (51)  and  yields,  using 
the  expression  Cl  for  D.(p), 


dJ(Q) 
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-i£_ 

(2  tt)3 
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'Pi 


(2  irY 


D.(Q  -  p) D.(p) 


°i(^l  "  Pi>D.(pi)27r6 


73 


PREVIOUS  PAGE 
IS  BLANK 


9 


74 


where  we  have  used  the  two  projection  delta  functions  to  integrate  over  q  and  p  . 
The  argument  of  &  (p^  -  q^)  may  be  expressed  using  Eq.  (62)  and  letting 
s  =  sin  ($f  -  c  =  cos  (${  -  <p.)  as 


r-*  -  n  )2 

q,  -P,  .  =  [-q9s  +  q  c  -  p  ]' 

L1  1 -I  along  p.  )  1  1  1 


_*  q  2  2  2  2 

q1  -  px  =  qi +  q2  ■  '  q2s>  ■  2P2*qis  +  q2c)  +  Pi 


This  yields 


nj  _  1  ,rLj.Li  ?rLiLf  f"  ttL2n4-  L(| 


4  (n  -  j) 


x//dPldqi  e^i-Ki/^/4fa-J»e-^i-KfA2/4J 

^  -Up1-K.i)tanx.i-Kcosxi-(q1-Kf^)tanxf+Kcosx.f]2L2/4 


X  /  /  dp„  dq„  e 


-(p2-Ki2)2L2/4(n-j)  -(q2-Kf2)L2/4j 


,"[q2l+Pl"2p1q1C]Li/4  y  -[q2+P2-2p2(q1S+q2c)+2Plq2s]L2/4 


Performing  the  integration  over  p„  and  q,  gives 


2  +2 

2,  2  , 


f  I"  -P2L2/4(n-j)  -q2  L^/4j 

Jdp2./dq2e  e  2  1 


X  e 


'  [^2+p2"2p2(qiS+q2C>+2plq2S]  Li/4 


4tt 

L2 


j(n  -  j) 


;_2, 


i  |_  1  +  n  +  js  (n  -  j) 


1/2  +p2s2L2/4(l+l/j) 


X  e 


[ -q^s+p jS- c/ (1+l/j)]  2/(?2 


(D6) 


where 


<PQ  =  [(n+  1)  +  js2 (n  -  j)]/t(n  -  j)(j  +  l)]  . 

Performing  the  Pj,  qj  integration  tl  en  yields 

i 


nj  r  L.  L.  L 

•  _  V£  _ l  f  || 


V^n  +  1  +  j(n  -  j)s2  ^  a2li2  +  y2 
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X  e 


to  j+fl  1'Y/2/32  ]  2/4[a2+72  /402  ]  02/4 0 


2  r 

e  e 


(D7) 


(D8) 


where 


2 

s  x. 


» 1  =  J-  \  -X.  +  cy  +  p2  (x  t2  -  y  t  t  )  +  - —  + 

1  “I  1  1  00  Joos 


2  2 

CSX, 


(1  +  l/j)  ^‘(l  +  l/j) 


cs  y, 


^o(1  +  l/j)J 

<7  =  i  (<>rn)+  (r)2to  +  ])  -ttmtu 


•>  ^ 
c“  s“ 


^  ;(i  +  i/j)“ 
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Appendix  E 

Multiscattering  Term  for  Gaussian  Correlation  Function  for  j  =  0,  n  -  j  >  0 

The  general  term  of  the  sum  of  Eq.  (54)  with  n'  =  (n  -  j)  >  0,  j  =  0  or  with  an 
interchange  of  subscripts  i  — —  f  for  (n  -  j)  =  0,  j  >  0)  can  be  written  for  the  Gaus¬ 
sian  correlation  function  as 


=  J  eiAI^'  Ar  ®  ^  (Ar)  D?'(Ar)  dAl- 


gip-  A~r  £  {?) 

J  (2*r 


dAr  e 


*  j 

J  (2  *r 

1  3/d  P 

/  dq 

(2»)  J 

l 

X  6  (K  - 

+  p  -• 

dq  *0<p)  D"  (qx)  •  6(q  •  K.) 
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Letting  c  =  cos  ( <p ^  -  <pj  anti  s  =  sin  (<?£  -  p.)  this  may  he  written  as 


*  /  “ 
J  C 


(2»)3K. 


*Li  Lll  vLlh 


1  1 


l2irK.  I 
n  iz 


-(Kfi  •  c-c,1)2l|/4  -(Kfis-q2)2L2/4 


-(q1-K.i)“L“/4n'  -(q2)2L2/4n' 


(  q  .  K.  -  K2  \2  „ 

i 


Carrying  out  the  indicated  integration  this  becomes 


2  2 

,n-  »/»l2lm 

'  =  ‘(n*  +  TT  e  L 


1  9  9  -“S' Ik  1-y ,c)-(x  -y  )n2tan(x.)] 

-KV’  *e\  1  1 

X  e  e 


where 


X1  =  KilLi  *  yl  =  KfiLl 


xa  «  K  cos  (X.)  Lx  ,  Kcos(  w  , 

a  IX 

r  9  2  2  ll/2 

Lt  =  L“  +  Ljj  tan  (Xj) 


q.  =  L./ L  , 

l  l  1 


and  where 


tP .  =  ll  +  n./n'  +  q2  tan2  (\.)l  1/?2 


til 


!-iep>e>ent<itiori  .»f  K'llmojoruH  «;orrui.r.iO) 


.juiv. 


In  order  to  appi  uj.'tj  ) 

that  it  cati  be  repra  sent oc!  in  a  form 


;  s  a  ,-mn  -  '  .  . ; . .  1 1  •.!>]'  l  i.i  ns. -a  ;■ ;  t.-'  ,  v 
similar  U*  O .  ).  where  from  L«, 


< 
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with 


(F10) 


CF 11) 


The  set  of  M  simultaneous  linear  equations  given  by  Hq.  (F7)  can  then  be  solved 
for  C<n). 

I 
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